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ABSTRACT 



We present the results of an observational photo-polarimetry campaign of the blazar PG 1553 + 1 13 at optical wavelengths. The blazar 
was recently detected at very high energies (> 100 GeV) by the H.E.S.S and MAGIC y-ray Cherenkov telescopes. 
Our high-temporal resolution data show significant variations in the linear polarization percentage and position angle at inter-night 
time-scales, while at shorter (intra-night) time-scales both parameters varied less significantly, if at all. Changes in the polarization 
angle seem to be common in y-ray emitting blazars. Simultaneous differential photometry (through the B and R bands) shows no 
significant variability in the total optical flux. We provide B and R magnitudes, along with a finding chart, for a set of field stars 
suitable for differential photometry. 
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1. Introduction 

The active galactic nucleus (AGN) PG 1553+113 was first 
detected in the P alomar-Green survey of UV bright objects 
dGreen etaD[l986). It is classified as a BLLac object based on 
its featureless spe ctrum dMiller & Green|[l983l) and significant 
optical variability (Mill er et al.lll988l ). The object has been well 
studied from radio to VHE y-ray (> 100 GeV), including sev- 
eral simultaneous multi-wavelength observatio n campaigns (see 
lOsterman eTaill200fjt iMankuzhivil et al.ll2009l) . 

At radio frequencies between 4.8 GHz and 14.5 GHz, the 
source was found to be variable on ti mescales of months 
dPerlman et al]2005l:IOsterman et al.l2006l) . VLBA observations 
have shown a jet extending at le ast 20 pc towards the northeast 
of the object dRector et al.ll2003l) . 

In the X-ray band, PG 1553+1 13 has be en detected at differ - 
ent flux lev els by Einstein ROS AT, ASC A (iDonato et al.ll200lb 



BeppoSAX (IDonato et alj|2005l) RXT E (lOsterman et al 



2006) 



XMM -Newton dPerlman et al.1 120051), Swift dtramacer e et all 
120071) . and Suzaku dReimer et al.ll2008l) satellites. A list contain- 
ing the flux and spectral para meters from all X-ra y observations 
of the source can be found in lAbdo et al.l d201 0b). 

At very high energy (VHE) the object was first detected with 
the HESS stereoscopic a rray of imaging atmosph eric-Cherenkov 
telescopes during 2005 ( Ah aronian et al.ll2.Q06b . This detection 
was l ater confirmed thro ugh MAGIC observations in 2005 and 
2006 (lAibert et al.ll2007l) . Both VHE spectra are unusually soft, 



* Based on observations collected at the Centra Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck Institut fur Astronomie and the Instituto de Astrofisica de 
Andalucfa (CSIC). 



although these results should be taken with caution because the 
errors are large. PG 1553 + 113 was also de tected at GeV y-ra y 
by the Fermi Gamma-ray Space Telescope (Abdo et al . 2010a). 

In this frequency range the source presented one of the 
hardes t spectra among the 106 AGNs listed in Abd o et al.l 
(2009). Combining archival radio, optical, X-ray, and VHE y- 
ray data, these authors modelled its spectral energy distribu- 
tion and demonstrated that a simple, one-zone synchrotron self- 
Compton (SSC) model provides a reasonably good fit of the ob- 
servational results. Unfortunately, the redshift of the source is 
still unknown. At present no emission or absorption lines have 
been measured despite several observational campaigns with op- 
tical instru ments. However, an uppe r limit of z < 0.74 was de- 
termined by lAharonian et al.l (120061) from the photon spectrum 
obtained by H ESS and later con firmed from MAGIC spectral 
measurements (Albert et al. 2 0071) . 

In this paper, we report optical polarimetric observations of 
the y-ray blazar PG 1553 + 113. Our high-temporal resolution 
data show significant variations in the position angle at intra- 
night as well as at inter-night time-scales. Simultaneous differ- 
ential photometry (at the B and R bands) shows no significant 
variability in the total optical flux. We describe the observations 
and data analysis in Sect. [2] Then, we present the obtained re- 
sults in Sect. [3] and we discuss the origin of the polarization 
angle variations and state our conclusions in Sect. |4] 



2. Observations and data reduction 

We observed PG 1553 + 1 13 during six nights in April 2009 with 
the Calar Alto Faint Object Spectrograph (CAFOS) in imaging 
polarimetry mode, at the CAHA 2.2 m. telescope, Calar Alto, 
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Spain. The detector was a 2048 x 2048 px blue sensitive SITe 
CCD, with 24/itn pixels (the resulting scale is 0.53 arcsec px). 
In this mode, CAFOS is equipped with a Wollaston prism plus 
a rotatable half-wave plate, which produce two orthogonally po- 
larized images -i.e., corresponding to the ordinary (O) and the 
extraordinary (E) beams- on the focal plane, with an effective 
beam separation of ~ 20arcsec. A mask with alternating blind 
and clear stripes is placed before the detector to avoid image 
overlapping. Four frames, each with a different position of the 
half-wave plate (0, 22.5, 45, and 67.5 deg) are needed to obtain 
the normalised Stokes pa rameters (U, Q) for linear polariza- 
tion (for more details, see ICellone et al.l 2007bl and references 
therein). 

We observed through B and R Johnson-Cousins filters. The 
images range in exposure time from 150 to 240 s and from 50 
to 120 s for the B and R filters, respectively. Photometric data 
are obtained by adding up the fluxes of the O and E images. 
In this way were able to follow simultaneously the polarimetric 
and the photometric behaviour of the source. Standard stars from 
iTurnshek et al.l (ll99Qh were observed to determine the zero point 
for the polarization angle and the instrumental polarization. The 
absolute calibration of the pho tometry was pe rformed by observ- 
ing several standard stars from Landolt (1992) over a wide range 
of air masses during photometric nights. For the same purpose, 
several images of the PG 1553+1 13 field were taken without the 
polarizing unit and through the B and R filters. The images were 
debiased and flat-fielded in the standard way using the IRAlQ 
reduction package. Regrettably, the whole first night and part of 
the last two nights were cloudy or with poor seeing conditions. 
We thus present results for the five clear/useful nights from April 
21 to April 25, 2009. 

Instrumental magnitudes corresponding to the O and E im- 
ages were obtained by means of aperture photometry using the 
IRAF task apphot. The same process was performed for nine 
field stars suitably placed with respect to the mask. We used 
these field stars to estimate the foreground polarization, while 
two of them were selected to perform the differential photome- 
try (see lCellone et al. 2007a for prescriptions on the appropriate 
choice of stars for differential photometry). In all the cases, a 
3 arcsec aperture was used. 

We obtained the linear polarization (P) and the position an- 
gle (6) through each filter for PG 1553 + 113 and the field stars 
along each night. The data were corrected for instrumental polar- 
ization using the results for the unpolarized standard stars. The 
foreground polarization (formally, a lower limit) was estimated 
from a suitable field star (#9, see Fig.[TJ, lying at 156 arcsec pro- 
jected distance from PG 1553+113, resulting in P fg = 0.24%. 
This value agrees with the low Galactic extinction in the direc- 
tion to t he blazar (E(B-V) = 0.052 mag. => P < 9 x E(B - V) = 
0.468%, Schle gel et al.l fl998). As a check, we measured the de- 
gree of polarization of several other field stars. The mean value 
in both filters, B and R, is < 1%, which confirms a low fore- 
ground polarization. The position angle was transformed to the 
standard system using the highly polarized standard stars. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Finding chart for comparison stars in the field of 
PG 1553+113. The image is 9 arcmin on a side, with north up 
and east to the left. The blazar is bracketed by two vertical bars. 

Table 1. Standard magnitudes for field stars. 



star 


B 


eB 


R 


eR 


# 


mag 


mag 


mag 


mag 


1 


14.52 


0.02 


13.45 


0.02 


2 


15.62 


0.02 


14.16 


0.02 


3 


16.45 


0.03 


15.26 


0.02 


4 


14.58 


0.02 


13.59 


0.02 


5 


13.95 


0.03 


12.74 


0.03 


6 


16.52 


0.03 


15.18 


0.02 


7 


14.66 


0.02 


13.58 


0.02 


8 


13.53 


0.02 


12.50 


0.02 


9 


15.92 


0.02 


14.98 


0.02 


10 


16.62 


0.02 


15.33 


0.02 


11 


15.75 


0.02 


13.79 


0.02 


12 


14.71 


0.02 


13.76 


0.02 


13 


15.20 


0.02 


14.00 


0.02 


14 


16.68 


0.03 


15.38 


0.02 


15 


16.84 


0.03 


15.77 


0.03 


16 


16.80 


0.03 


15.95 


0.03 


17 


16.62 


0.03 


15.58 


0.03 



3. Results 

3.1. Photometry 

The differential light-curves were obtained in the usual way, us- 
ing a non-variable star in the field as comparison, while another 
(also non-variable) star was used to construct a second differen- 
tial light-curve against the comparison star, to be used for control 
purposes. Fig.[T]is a finding chart for non-variable stars Q suitable 
for differential photometry in the fi eld of PG 155 3+113. Their 
magnitudes in the standard system dLandoltll 19921) are given in 
Table Q] Throughout this study we used stars 12 and 11 as com- 
parison and control, respectively. 



2 The stability of these stars was verified with the 2007 data, see be- 
low. 
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A statistical analysis of these light-curves reveals that this 
VHE blazar showed no significant variability in B or in R. This is 
true considering each individual night (microvariability) as well 
as taking the five nights together (inter- night variabili t y). The 



as taking tne five ntgnts together (inter-nignt variability 
statistical analysis was performed as in ICellone et al. (C 



(2007. 



and applying the corrective factor F (see iHowell et al] 1198 
which weights the statistics taking into account the differences 
in magnitude among the target and the comparison and control 
stars. The values for T were ~ 0.75 for B and ~ 1.12 for /fl 

The dispersions of the curves ranged from 0.004 to 
0.013 mag, and the mean apparent magnitudes of PG 1553+113 
for the whole campaign were (B) = 14.90 + 0.003 mag and 
(R) = 14.18 + 0.003 mag. Differential photometry data taken 
with the 2.15-m "Jorge Sahade" telescope, at CASLEO, San 
Juan, Argentina, during April 2007 seem to point in the same 
direction: the source was non-variable at optical wavelengths. 
Because of the mask used for CAHA observations, different 
comparison and control stars were used for both data sets, but we 
were able to put all results on the same magnitude scale because 
they had several field stars in common. Although the blazar's 
flux remained steady during each of both periods, the object was 
about 0.1 mag fainter in the /?-band during 2009 compared with 
2007. The results on PG 1553+1 13 from the April 2007 observ- 
ing run, along with those for other TeV y-ray emitting blazars, 
will be presented elsewhere (Andruchow et al., in preparation). 

Both observations were perf orm ed after the high-en ergy de- 
tections dAharonian et al.ll2006l and lAlbert et al.ll2007l) . On the 
other hand, the observations reported here were almost simul- 
taneous with the final part of the first ~ 200 days of the 
Fermi Gamma-ray Space Telescope science operations (ended 
on February 22, 2009), when PG 1553 + 113 was stud ied in de- 
tail at the GeV gamma-ray regime (Ab do et al.l20 10a). The GeV 
emission state of PG 1553+113 did not show any significant 
changes; this can be taken as an indication that it was in a low- 
energy state throughout the observations. 

The photometric spectral index taking the five nights to- 
gether and using the mea n magnitudes (corrected for Galactic 
extinction, ISchlegel et al.l 119981) was a BR = -0.71 + 0.01. 
iFalomo et al.l dl994) reported a mean value for the spectral index 
at visual wavelengths of (a) = -0.96, obtained from continuum 
fitting to the optical spectrum. Part of the difference between the 
two values could arise because our a was calculated using pho- 
tometric observations, which could be affected by atmospherics 
effects not completely removed (e.g., telluric lines). During the 
period of five years in the early '90s, the authors reported that 
the maximum variation in a was of about 0.24. 



3.2. Polarimetry 

Figure |2]presents the curves for the linear polarization and posi- 
tion angle against time along the five nights. 

Following a x 2 criterion, we performed the statistical anal- 
ysis for each night, in each filter, and also for the whole cam- 
paign. According to this criterion, a source can be classified as 
variable in an observing session for the observable S if the prob- 
ability of exceeding the value X 2 by chance is < 0.1 %, and it 
is classified as non-variable if the probability is > 0.5 %. If no 
definite decision could be reached with that criterion (any value 
between 0.1% and 0.5%), we call this a dubious behaviour. If 
the errors are random, X 2 should be distributed as x 1 with n - 1 
degrees of freedom, where n is the number of points in the distri- 



3 T = 1 for the ideal case when target and both stars all have the same 
magnitude. 



bution. In Table [2] we show the values of the variability param- 
eters for the linear polarization percentage (upper section) and 
for the position angle (lower section). Column 1 gives the obser- 
vation date, Col. 2 gives the corresponding filter, Col. 3 shows 
the number of points for each night, Col. 4 gives the mean po- 
larization (position angle) for the observing night, Col. 5 shows 
the rms o~p (o~g), Col. 6 gives the value for AP (A9) (the differ- 
ence between the maximum and the minimum values during the 
period reported), Col. 7 is At (between the maximum and the 
mi nimum), and Col. 8 show s the value of^- 2 . 

lAndruchow et al] (120081) present a method to estimate the 
degree of polarization of the nucleus, corrected for the de- 
polarizing effect from the host galaxy light. This requires a de- 
tailed knowledge of the host's absolute magnitude and effective 
radius, but this information is not available forPG 1553 + 1 13 (its 
host remains undete cted). However, fr om the upper limits on the 
magnitudes given in lUrrv et al.l (120001) . we do not expect that the 
host galaxy of PG 1553+1 13 would significantly reduce the de- 
gree of polarization of the active nucleus, nor would it introduce 
any spurious variation. 

4. Discussion 

We found that the flux in both filters, B and R, remained steady 
during the observing run reported here. On the other hand, for 
the degree of polarization and the position angle in the same 
period and in the same filters, the source showed small varia- 
tions. Comparing our results with previous ones, we remark that 
lOsterman et al.l d2006l) found in their multi wavelength campaign 
(April-May 2003) that PG 1553+113 was brighter in both op- 
tical bands by about 0.5 mag without signs of variation. Those 
obse rvations were made durin g an X-ray flare. In the same pe- 
riod, [Andruchowetal] (|2005) followed the source in linear op- 
tical polarization (V filter) during two nights with a good sam- 
pling in both nights. The degree of polarization and the position 
angle then showed microvariations with a most likely flickering 
behaviour. For the values reported here (see Table |2}, we note 
that the degree of polarization was higher during 2003, as was 
the position angle. The polarization vector seems to have rotated, 
because the position angle varied by ~ 20° between both epochs, 
but it wa s still roughly ortho gonal to the radio-jet direction re- 
ported by iRector et all d2003l) . 

This is particularly interesting given the recent detection of 
rapid optical polarization cha nges correlated wi th y-ray vari- 
ability in the blazar 3C279 dAbdo et al] 1201 Od) . which have 
been interpreted as the result of geometrical effects (see also 
lAndruchow et~aT1l2003l) . 

We also have preliminary results f rom photometric and po- 
larimetric data taken in April 2007 (lAndruchow et al.l [2007) 
trough R and V filters. These show that PG 1553+1 13 presented 
no flux variations during the four nights during which it was fol- 
lowed, although the blazar was then slightly brighter ( ~ 0. 1 mag) 
than in 2009. Regarding the optical linear polarization, we ob- 
tained one point per night, with no significant variations. The 
mean value was (P) = 2.77%. 

During 2005 and 2006 this sou rce was detected at v ery 
high energies by H ESS and MAGIC dAharonian et al]|2006l and 
lAlbert et aD 120071 respectively). This very high emission was 
detected in between the optical observations that we discused 
above. 

The data presented here can be useful to complete the multi- 
wavelength view of this peculiar object. In the near future, using 
deep exposures with the Gemini North 8 m telescope, we will 
attempt to determine the redshift of the host galaxy. 
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Table 2. Statistical results for the linear polarization. 



Date 


Filter 


n 


(P) 


(Tp 


AP 


At 


x 2 


(2009) 






% 


% 


% 


hr 




21 Apr 


B 


7 


1.92 


0.27 


0.80 


3.56 


13.96 (4) 


21 Apr 


R 


7 


1.95 


0.15 


0.40 


3.05 


8.78 (4) 


22 Apr 


B 


10 


1.80 


0.14 


0.41 


4.94 


8.77(1) 


22 Apr 


R 


10 


1.82 


0.16 


0.56 


4.47 


10.30 (3) 


23 Apr 


g 


13 


2.05 


0.16 


0.53 


2.72 


11.87 (8) 


23 Apr 


R 


14 


2.00 


0.11 


0.42 


4.37 


6.99 (1) 


24 Apr 


B 


2 


2.36 


0.36 


0.51 


0.60 


1.91 (8) 


24 Apr 


R 


2 


2.30 


0.26 


0.37 


0.57 


1.66 (6) 


25 Apr 


B 


2 


3.06 


0.38 


0.53 


0.62 


1.07 (4) 


25 Apr 


R 


3 


2.86 


0.28 


0.50 


1.35 


1.38(4) 


All 


B 


34 


2.03 


0.35 


1.95 


88.73 


80.35 (1) 


All 


R 


36 


2.03 


0.32 


1.75 


69.44 


77.61 (8) 




Date 


Filter 




w 


ere 


AO 


At 


x 1 


(2009) 












hr 




21 Apr 


B 


7 


120.24 


2.06 


5.11 


2.11 


7.06 (4) 


21 Apr 


R 


7 


116.82 


1.82 


5.76 


3.00 


7.52 (7) 


22 Apr 


B 


10 


131.75 


1.74 


5.09 


4.97 


5.86 (4) 


22 Apr 


R 


10 


129.27 


1.65 


4.98 


5.54 


6.99 (3) 


23 Apr 


B 


13 


129.86 


1.58 


5.44 


5.77 


5.76 (5) 


23 Apr 


R 


14 


129.75 


1.80 


5.38 


1.25 


10.19(2) 


24 Apr 


B 


2 


128.77 


0.20 


0.28 


0.60 


0.00 (4) 


24 Apr 


R 


2 


132.89 


0.54 


0.76 


0.57 


0.04 (8) 


25 Apr 


B 


2 


124.19 


2.51 


3.55 


0.62 


0.53 (9) 


25 Apr 


B 


3 


126.98 


0.76 


1.37 


0.63 


0.22 (2) 


All 


B 


34 


128.04 


4.69 


15.88 


26.07 


217.18 (4) 


All 


R 


36 


127.05 


5.45 


19.34 


69.75 


361.33 (1) 
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ter (left), and through the R filter (right). 



Osterman, M. A., Miller, H. R., Campbell, A. M., et al. 2006, AJ, 132, 873 
Perlman, E. S., Madejski, G, Georganopoulos, M., et al. 2005, ApJ, 625, 727 
Rector, T. A., Gabuzda, D. C, & Stocke, J. T. 2003, AJ, 125, 1060 
Reimer, A., Costamante, L., Madejski, G, Reimer, O., & Domer, D. 2008, ApJ, 
682, 775 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Tramacere, A., Giommi, P., Massaro, E., et al. 2007, A&A, 467, 501 
Turnshek, D. A., Bohlin, R. C, Williamson, R. L., et al. 1990, AJ, 99, 1243 
Urry, C. M., Scarpa, R., O'Dowd, M., et al. 2000, ApJ, 532, 816 



4 



